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■ Abstract 

o, 

I The atomic cascade in fi p and vr p atoms has been studied with the improved 

version of the extended cascade model in which new quantum mechanical calculations 
of the differential and integral cross sections of the elastic scattering, Stark transitions 
Q i and Coulomb de-excitation have been included for the principal quantum number 

values n ^ 8 and the relative energies E ^ 0.01 eV. The X-ray yields and kinetic 
energy distributions are compared with the experimental data. 

1 Introduction 

o 

^ ' The exotic hydrogen-like atoms are formed in highly excited states, when a heavy negative 
^ ■ particle vr", K~, etc.) slows down and captures in hydrogen media due to inelastic 
collisions. The formation is followed by the atomic cascade in which the initial distributions 
, in the quantum numbers (n, I) and kinetic energy of the exotic atom change due to various 
O ! processes: radiative transitions, external Auger effect, elastic scattering. Stark transitions. 
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Coulomb de-excitation, weak decay and nuclear absorption (in case of hadronic atoms). A 
' good understanding of the kinetics of atomic cascade in hydrogen-like exotic atoms is very 
^ ■ important both for the planning and interpretation of the experiments (see for a review 
^ . [1]). The first theoretical study of the atomic cascade was performed more than forty five 
[ years ago by Leon and Bethe 0. In this paper, and later in more refined models [3l H], 
the rates of the collisional processes were calculated in the semiclassical approximation at 
a fixed kinetic energy (~ 1 eV) and used for the simulation of the atomic cascade. A 
number of experiments in which the energy distributions of fi~p, fi~d and n'p atoms were 
measured using various time-of-flight methods [H [6l [3 El E] showed that the kinetic energy 
changes during the cascade. In particular, the existence of high energy components has been 
established in pionic hydrogen in the neutron time-of-flight experiment [7J and in muonic 
hydrogen in diffusion experiments [Bl E]- This experiments require a more sophisticated 
approach taking the evolution of the kinetic energy distribution into account during the 
whole atomic cascade. Beginning with Markushin [10] recent cascade models take into 
account the acceleration and deceleration processes. The extended standard cascade model 
(ESCM) [11] introduces a number of improvements compared to the earlier models: for 
example, the scattering from molecular hydrogen at high n is calculated as opposed to the 
phenomenological treatment in other cascade models. 

In the recent papers [T2t [13] , the dynamics of the collisions of the excited exotic atoms 
with hydrogen ones has been studied in the framework of the close-coupling approach. The 
elastic scattering. Stark transitions and Coulomb de-excitation were described by a unified 
manner in this approach. The differential and integral cross-sections were calculated in a 
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wide range of the principle quantum number values and relative energies for muonic, pionic 
and antiprotonic hydrogen atoms. The vacuum polarization and strong interaction shifts 
(for hadronic atoms) were taken into account (see also this issue). 

In the present paper the first theoretical results of the X-ray yields and kinetic energy 
distributions for fi~p and 7r~p atoms calculated in the new version of the ESCM are intro- 
duced. The new results for the coUisional processes ([I2lll3j, see also this issue), especially 
for Coulomb de-excitation process, are used in the present ESCM and lead to an essential 
improvement over the previous cascade calculations. 

2 The extended standard cascade model 

The cascade in exotic atoms (as in the previous ESCM) is divided into the classical domain 
for high n and the quantum mechanical domain for low n. In the classical domain (for 
muonic and pionic hydrogen we use n > 9) the classical-trajectory Monte Carlo calculations 
have been included in the cascade model (see for details [11]) with the molecular structure 
of the target taken into account for the description of the elastic scattering. Stark mixing 
and Coulomb de-excitation processes. The external Auger effect was calculated in the 
semiclassical approximation through the whole cascade. In the upper stage of the cascade 
the distributions in quantum numbers and kinetic energy are calculated. Further, these 
distributions were used as input data for the next stage of the atomic cascade. In the 
quantum mechanical domain (at n < 8 ) the differential and integral cross sections for all 
nl I— > n'l' transitions (besides the Auger and radiative de-excitations) were calculated in the 
fully quantum mechanical close-coupling approach (see [121 [13] and references therein) for 
the initial states with n < 8 and relative energies E > 0.01 eV. The corresponding arrays 
of the cross sections are included in the cascade code. Thus the present cascade code does 
not employ any fitting parameters and the kinetics of the atomic cascade is treated more 
accurately. 

The initial conditions for the cascade calculations are defined by the initial distributions in 
the quantum numbers n and / and the laboratory kinetic energy of the exotic atoms. In 
the simplest picture, the exotic particle is captured by the proton of the hydrogen atom 
in a state with ~ \/M (M is a reduced mass of exotic atom in atomic units). More 
elaborate approaches [HI |15] take the molecular effects of the target into account and 
predict distribution in the initial rii with the peak shifted towards lower values and non- 
statistical Zj-distribution. In the present cascade calculations the influence of the initial 
parameters are also studied. In order to obtain good statistics the cascade calculations 
were performed with 10^ events. 

3 Results 

3.1 Muonic hydrogen 

Muonic hydrogen atom is the simplest of exotic atoms and its experimental and theoretical 
studies are the best probe for the investigation of the various processes in exotic atom from 
hydrogen atom or molecule collisions. The present cascade calculations in muonic hydro- 
gen have been done in the density range (10^^ — 1) (in units of Liquid Hydrogen Density, 
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LHD = 4.25 • 10'^'^ atoms / cm^) using various initial conditions. The results are compared 
with the experimental data [T^ IT7t [T5] for relative K X-ray yields. The energy distribu- 
tions of muonic hydrogen and Doppler broadening of the X-ray lines are also introduced. 




Figure 1: The density dependence of the absolute K X-ray yields in muonic hydrogen. 

In contrast to the practically linear dependence of the absolute Ktot yield on the density 
in the wide range from 10^'' LHD to ~ 10^^ LHD the individual Ki yields [i = a, (3, etc.) 
have a more complicated density dependence due to the competition between the various 
coUisional and non-coUisional processes during the cascade. The atomic cascade at low 
target densities is mainly dominated by the radiative de-excitation proceeding through the 
circular states that result in 2p Is transition whereas the transitions np Is are much 
weaker. We see the corresponding picture in Figure 1. At very low density the kinetic 
energy distribution of the muonic atom must conserve it initial kinetic energy distribution. 
While the density is increased the role of the collisional processes enhances leading to both 
the decrease of Ka yield and simultaneously the increase of the other Ki yields {i ^ a), 
conserving the total absolute yield practically unchangeable. With increasing density the 
Coulomb de-excitation becomes also more efficient (especially at low energies) and leads 
to the acceleration of the exotic atom and higher populations of the np states with n > 2 
due to the Stark transitions. At the target densities more than ~ 10~^ -r- 10~^ LHD the 
collisional de-excitation n —>■ n'{n' < n) begins to dominate the radiative transitions (ex- 
ternal Auger effect and Coulomb transitions) resulting in the decrease of the K>p yields. 
Since the collisional de-excitations to the Is state of the muonic hydrogen are strongly 
suppressed (more or about 2 keV energy is released) we observe the enhancement of the 
absolute Ka yield due to the collisional transitions to 2p state and the Stark 2s 2p tran- 
sition at the energy above the threshold. Thus the swift suppression of the absolute Ka 
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yield at the densities from 10"'' LHD to ~ 10"^ LHD and increase at more higher densities 
can be explained by the competition of the different processes during the atomic cascade 
in which the evolution of the kinetic energy is taken into account through the whole cascade. 

The discussed above regularities in the absolute yields can be justified or vice versa in 
comparison of the calculated relative X-ray yields for the iJ,~p with the experimental data 
[T6l \T7\ [T8] introduced in Figure 2. 




N/LHD N/LHD 



Figure 2: The density dependence of the relative X-ray yields in muonic hydrogen: Ka 
on the left; Kf^, and yields (black, red and green, respectively) on the right. The 
experimental data are from [161 [13 IH]- 

The agreement between theoretical results and experimental data is very good practi- 
cally for all densities under consideration. The observable disagreements can be a subject 
for further both theoretical and experimental studies. In our opinion, there are also a num- 
ber of problems which demand more elaborate experimental studies: a proper separation 
of the different k lines, inclusion a Doppler broadening effect in the analysis of the experi- 
mental data and so on. It would be very important to check the present results directly by 
measuring the absolute X-ray yields at some values of the target density. 

The kinetic energy distribution of the exotic atom changes during the cascade and is a 
more refined probe of the theoretical approaches to the description of the cascade processes. 
Using the cascade model described above we calculated the kinetic energy distribution 
at the moment of the i^Q,, Kp, etc. radiative transitions and the corresponding Doppler 
broadening of the Is line due to the kinetic energy distribution of the exotic atom at the 
instant of the np — > Is radiative transitions (ra = 2 -^ 4) at the target pressure 10 bar 
(~ 10"^ LHD). The results are shown in Figure 3. 

As it is seen from Figure 3, the kinetic energy distribution has distinctive high energy 
structures arising from different Coulomb transitions with An > 1 preceding the radiative 
de-excitation. The complicated shapes of these structures can be simply explained by the 
interplay of three factors: the kinetic energy distribution of exotic atom before Coulomb 
transition, the anisotropy of the angular distribution in the Coulomb de-excitation pro- 
cess and, finally, the deceleration due to the elastic scattering and Stark transitions after 
Coulomb de-excitation but before the radiative de-excitation. 

The energy distribution at low n can be determined by measuring the Doppler broadening 
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Figure 3: On the left: the kinetic-energy distribution of jj,~p at the instant of the radiative 
np — >■ Is transitions in gaseous hydrogen at pressure 10 bar. On the right: the i^T-hne 
shapes after Doppler broadening due to the kinetic-energy distribution. 



of K X-ray hues. As demonstrated in Figure 3, the high-energy components lead to the 
significant Doppler broadening, which especially pronounced for line. The resulting 
shapes and widths for Ka,Kp and lines at target pressure 10 bar are shown in Fig. 3 
(on the right panel). A significant spreading up to ±1.25 eV for line was found. 

The essential acceleration during cascade is illustrated in Figure 4. Here the calculated 
density dependence of the mean kinetic energy at the instant of the radiative np Is 
de-excitation is shown. At the densities from 10^'' LHD to ~ 2 - 10~^ LHD the mean energy 
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Figure 4: The density dependence of the mean kinetic energy of the fi p at the instant of 
the different radiative np ^ Is transitions. 
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slowly increases for all lines and changes correspondingly from 4 eV to 8 eV. At more higher 
densities the mean energy of the Kp line reaches the maximum value ~ 20 eV in the density 
interval 10~^ — 10~^ LHD. The most pronounced increase of the mean energy was obtained 
for Kd line with the maximum more than 100 eV at the density ~ 8 ■ 10~^ LHD. 

3.2 Pionic hydrogen 

The present study of the atomic cascade in pionic hydrogen was focused on the experimental 
data for the absolute X ray yields [I9] and the neutron time-of-flight spectrum [7J obtained 
at LHD. The results of our cascade calculations for the absolute X-ray yields as a function 
of the target density are shown in Figure 5 in comparison with the experimental data [19] . 
The nuclear absorption in the hadronic atoms results to a strong suppression of the KX-iay 
yields, beginning from very low densities. As a whole the calculated yields are in a good 
agreement with the experimental data. In the calculation we introduced the additional 
parameter simulating the induced absorption from np states. The value of this parameter 
used in all present calculations for 2p-state was chosen equal to 0.045 meV. 




10"' lO"*^ 10"^ 10"'* 10'^ 10"^ 10"' lo' 

N/LHD 



Figure 5: The density dependence of the absolute X-ray yields in vr p atom. The experi- 
mental data are from p9]. 

The kinetic energy distribution of the pionic hydrogen at the moment of nuclear ab- 
sorption in liquid hydrogen is shown in Figure 6. The experimental data from PSI exper- 
iments [7] (Bad. and the present theoretical description are shown on the left side and 
right side of the Figure 7 respectively. The most promising features are that for the first 
time we can explain in more details both the whole kinetic energy distribution of the 7r~p 
atoms at the instant of charge-exchange reaction and fine structure of this distribution. In 
particular, according to our calculations the component around ~ 105 eV and ~ 209 eV 
are explained by the 5 — >■ 3 and 3 — > 2 Coulomb transitions, respectively. Besides, the 
relative yields of high energy components at the energies ~ 105 eV and ~ 209 eV are in 
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good agreement with the experimental data [7]. Thus in contrast to previous theoretical 
studies the prediction of the present cascade model are in good agreement with the neutron 
time-of-flight spectrum at LHD. 
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Figure 6: The kinetic energy distribution of TT~p at the instant of the nuclear absorption in 
liquid hydrogen: the experimental data [7] ( on the left) in comparison with the results of 
the present cascade calculations (on the right) 



In order to illustrate both the importance of the CD transitions with An > 1 and 
strength of the CD transitions with An = 1 the three variants of the cascade calculations 
were performed: full calculation where all CD transitions are taken into account, the cal- 
culation with only An = 1 CD transitions and, finally, the calculation with An = 1 CD 
transitions decreased by factor ken = 0.1. The results introduced in Figure 7 show that in 
contrast to the CC approach [12], |13] the previous theoretical description of the CD process 
in the semiclassical [20] and advanced adiabatic approaches (see [21] and references therein) 
are in disagreement with the experimental data [7]. 



0.02 
0.018 
0.016 
0.014 
0.012 

0.01 
0.008 
0.006 
0.004 
0.002 




i ' 1 


1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 
density LHD 


1 ' 1 


1 ' 1 

all An 

An=l 


1 ' 1 


- 
















T T 






T - 


1 


I.I 




1 , lill , 1 





20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 
E, eV 



0.02 
0.0 IS 
0.016 
0.014 
0.012 

0.01 
0.008 
0.006 
0.004 
0.002 




I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I ' I 



density LHD 



I J. T T 



T T 




all An 

<10 







.^lJ^^K.--f'l4^J_.. L ,-1 ^.JVll. -, .L ^J-J 



20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 
E, eV 



Figure 7: The kinetic energy distribution of iv^p at the instant of the nuclear absorption 
in liquid hydrogen for the different variants of the Coulomb transitions: full calculation - 
green line and only An = 1 CD transitions - blue line (on the left side); full calculation - 
green line and only An = 1 with the factor kcd = 0.1 - blue line (on the right side) 
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4 Conclusion 



The present version of the cascade model significantly improves the previous cascade cal- 
culations due to the self-consistent treatment of the elastic scattering, Stark transitions 
and Coulomb de-excitation. The inclusion in the cascade code the results of the detailed 
quantum mechanical calculations of the differential and integral cross sections for nl i— > n'l' 
collisional transitions allowed us to describe the kinetics of the atomic cascade without 
employing any fitting parameters and more accurately. 

The predictions of this cascade model are in a good agreement with the experimental data 
both for muonic and pionic hydrogen atoms. In particularly, the present cascade calcula- 
tions allow for the fist time to explain the observed kinetic energy distribution of 7r~p atoms 
at the instant of nuclear reaction. Our results explain the high energy components around 
~ 105 eV and ~ 209 eV (due to 5 ^ 3 and 3^2 Coulomb transitions) and lead to a very 
good agreement with the experimental weights of these components. 
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